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Abstract 

The heat capacity of an individual mesoscopic normal (nonsuperconducting) fermion 
system treated as a canonical ensemble of independent particles confined in a spa- 
tial area with fixed confinement parameters is studied in a wide range of particle 
numbers N < 10^ and temperatures which vary from values close to zero up to 
the Fermi energy sp- The temperature evolution of the heat capacity is naturally 
divided into four stages. On the first one the heat capacity exponentially increases 
with temperature and at a resonance temperature reaches either a local maximum 
or an irregularity in its growth. This temperature being measured can give informa- 
tion concerning the level spacings As in the immediate proximity ep- Calculations 
of Ae ( as the difference between the Fermi level and the level just above it) per- 
formed for diverse systems confined within nonsymmetric oscillator or rectangular 
wells show that nearly uniform distributions Ae is stretched up to ~ (2 'i)SF {^f 
is the averaged level spacing in the vicinity of the Fermi energy). On the second 
stage (T is more or of order of Ae) the progressive suppression of the level density 
oscillations takes place. During this stage the heat capacity oscillations v.s. parti- 
cle number can be distinctly observed. These oscillations give the particle number 
variation of the temperature averaged level density. The A^-oscillations discontinue 
for all 3D-systems at the temperature of order of ep/N^/^ and for 2Z)-systems at 
ep/N^^"^. The growth of the heat capacity on the third stage of the evolution is gov- 
erned by the T-linear law if the particle number is large enoug h (iV > 10^). In small 
systems (A^ < 10^) the heat capacity reveals marked deviations from this law imme- 
diately near the temperature suppressing the level density oscillations. It is found 
for rectangular potentials that the Sommerfeld factor a in the linear law (C = aT) 
for such small systems shows more complicated A^-dependence as compared with 
large systems where a is strictly proportional to N. For temperatures tending to 
ep the deviations from the T linearity are evident and at T > ep (the fourth stage 
of the evolution) any system attains to the classic Boltzmann-Maxwell limit irre- 
spective of the particle number. All our calculations are carried out by using the 
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canonical polynomial method. The results of approximate methods (in which the 
Fermi-Dirac function is used as the occupation probability and the chemical poten- 
tial A is temperature dependent to allow for the particle number conservation on the 
average) are shown for some exemplary cases which indicate that the approximation 
explicitly taking into account dX/dT is almost equivalent to the canonical method 
excepting a temperature region T < 6f- 

Key words: Mesoscopic systems; Canonical heat capacity; Low temperature 
resonances 

PACS: 65.80.+n; 73.22.Dj 



1 Introduction 

The first investigation of the mesoscopic fermion heat capacity belongs ap- 
parently to Frolich , Ref. [1] of 1937, who considered a system with a single 
electron spectrum similar to the spherical oscillator, i.e. with equal level spac- 
ing Ae and high level degenerations that was regarded as an approximation 
to the spectrum of a cubic rectangular potential. This model allowed Frolich 
to establish such distinctive feature of the low temperature heat capacity (C) 
as exponential growth with increasing temperature. He also found that at 
subsequent heating up to fc^T ~ 2Ae the specific heat gains the bulk value. 

The most important objection against Frolich's approach was given by Kubo 
(summary of the investigations of Kubo and his collaborators is given in 
Refs [2||3] ) who emphasized the irregularities of level spacings in mesoscopic 
systems which, as it seemed untill the eighties of the past century, could not 
have a defined and fixed shape and size. Therefore for describing an ensemble 
of such objects Kubo introduced into consideration the probability distribu- 
tion for the spacings between adjacent levels. In this case the thermodynamic 
quantities had to be calculated not with the spectrum of an individual system 
but with interlevel spacings which possess some probability of their appear- 
ance. As Kubo supposed the Poisson distribution in which very small spacings 
are most probable this immediately results in practical disappearance of the 
exponential part of the heat capacity and even at T/6 < 1 {6 is the averaged 
level spacing) C becomes linear in temperature. Application of other types 
of the level statistics remains rapid growth of C v.s. T though now it is not 
exponential increasing but proportional to either or [2j. 

However the interest in such description of mesoscopic systems has partly di- 
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minished after the discovery of the shell structure in metal clusters [1] and 
appearance of the possibility to fabricate mesoscopic systems such as artificial 
atoms or quantum dots with given geometry and controlled particle num- 
bers Electron-shell effects were also identified in experiments on diverse 
metal nanowires [71I81I9] . One more stimulus for the consideration of individual 
fermion systems came from experiments on trapped gases of fermions [10] . 

Therefore on this new stage in development of mesoscopic physics when variety 
of phenomena could be explained only if the level structure was explicitly 
taken into account the main interest has been concentrated on investigations 
of individual systems with a fixed particle number and geometry. 

As any of such systems can exchange energy with surroundings it has to be 
treated as a canonical ensemble. For the first time an analytical formalism 
(approximate in general but exact for very large particle number) was worked 
out by Denton, Miihlschlegel and Scalpino |TT] for ID-oscillator. Brack et 
al [T2] applied the canonical formalism (though not in its simplest version) to 
calculations of the electron heat capacity of alkali clusters with realistic single 
electron spectrum created by the temperature self-consistent potential. Par- 
allel with exact canonical calculations approximate methods have been devel- 
oped: partial projection, Refs. [T3fT^ . saddle-point approximation, Ref. [15j . 
integral-transformation method, Ref. [16] . In this chapter we apply our canon- 
ical polynomial formalism, Ref. [T7] , the general idea of which is similar to that 
of Denton et al however our method is applied to any energy level distribution 
and arbitrary particle numbers. 

Though among many experiments performed on mesoscopic systems to the 
present time there are only several investigations of the thermal behavior and 
in particular the heat capacity of normal and superconducting mesoscopic 
systems (see e.g. Ref. [TSj and references therein) we believe that such ex- 
periments will gain greater development owing to their practical applications. 
Therefore preliminary theoretical estimations of C for individual systems in 
simplified models that we give in this chapter can be a basis for arrangements 
of future experiments and further theoretical improvement. 

Our calculations of the fermion heat capacity are carried out in the indepen- 
dent particle model that can be directly applied only to trapped ultra-cold 
Fermi-gases. In the majority of the electron systems the strong interaction of 
each electron with ion surroundings and other electrons forms fermion quasi- 
particles with an effective mass different from the free electron mass and an 
energy spectrum similar to that for free electrons in a confinement potential. 
As shown in the Fermi liquid theory [19] the quasiparticle potential ( the ex- 
istence of which is directly established by experiments on the observation of 
the shell structure) can depend on the excitation energy and temperature that 
does not explicitly taken into account in our calculations - this is our first sim- 
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plification. The second simplification consists in neglecting the quasiparticle 
damping i.e. levels are considered to be perfectly sharp. However both these 
simplification can affect only high temperature values of C. Besides, we do 
not allow for the residual electron-phonon interaction i.e. such its part which 
in an averaged form does not enter into the electron mean field. As known [3] 
this interaction gives rise to renormalization of Sommerfeld's constant a for 
bulk metals at those temperatures that lead to the linear variation of C v.s. T 
(C = aT) therefore our estimations of C can be essentially corrected by tak- 
ing into account the electron-phonon interaction in this temperature region. 
However as shown below at these temperatures (when C is proportional T) 
the heat capacity practically does not bear information concerning the level 
structure near the Fermi energy and so we hope that the electron-phonon 
interaction being allowed for does not change drastically our conclusions. 

Two types of confinement are employed in our calculations. As a hard confine- 
ment we make use of the 2D and 3Z) rectangular potentials while the harmonic 
oscillator potentials are regarded as a soft confinement. The realistic case is 
something intermediate between these potentials. This can reveal itself again 
in the temperature region that provides the linear variation of C with in- 
creasing T. The finite depth of the potential can affect the high temperature 
variations of C which is certainly only of academic interest in the majority of 
cases. 

The content of this chapter includes 7 sections. In Sec. 2 the basic formulae 
for the calculation of C and interpretation of numerical results are given, here 
also the stages of the temperature evolution of C are outlined in consecutive 
order. Then in Sec. 3-6 we give the results corresponding to each stage. In Sec. 7 
we present conclusions. 



2 The heat capacity as an average of the level density 

Thermodynamic properties of fermion isolated mesoscopic systems that can 
be described by the independent or free particle model are determined by their 
single-particle energy spectra depending on the type of the fermion confine- 
ment and the size of the system. For such systems the heat capacity C can be 
calculated as an integral involving the temperature variation of the occupation 
numbers n(e), i.e. dn{e)/dT (the temperature T is measured in energy units), 
and the exact level density 
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Et > 0, dt > 1 are respectively the energy and degeneration of each single 
particle level t: 

C/ks = JpUe)e^de; (2) 



We shall consider C for a fixed particle number N that imposes on function 
n{e) the condition 

oo 

^ = J Pex{e)n{e)de. (3) 



Eq. ([3]) is fulfilled automatically if a fermion ensemble is treated as the canon- 
ical one {CE) and n{e) is the canonical occupation number. If n{e) is replaced 
by the Fermi-Dirac function f{e) 

/(e) = [1 + 6^(^-^)1 (3 = T-^ (4) 



then Eq. ([3]) defines the temperature dependent chemical potential A. 

The temperature independence of {dN/dp = 0), Eq. Q), in which also as 
in Eq. ([2]) function n{e) is replaced by f{e), results in the representation of C 
in the form 



oo 

c/ks = J pex{.s)ip{e)de, (5) 



<^(£) =[3\e-\- /5^)2e^(--^) [l + e^(^-^)l (6) 

op L J 



Eq. ([5]) allows the interpretation of the heat capacity as an average of the level 
density and the role of an averager is played by (y9(£:)-function. 

To stress the difference in the description with the temperature dependent 
A and a fixed N from the grand canonical ensemble description, in which 
the independent variable is A, we shall apply the term the equivalent grand 
canonical ensemble (EGCE) [2Q| . 

The simplest approach to the canonical description is the so called grand 
canonical ensemble (GCE) approximation. This approximation accepts f{e) 
as the occupation probabihty with temperature dependent A (defined by Eq. ([3])) 
but the temperature derivative of A {dX/dj3) is omitted in calculations of C 
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in contrast with the EGCE description. The GCE approximation gives the 
correct results at those temperatures when j3d\/dj3 is very small as compared 
with A, that e.g. leads to the correct linear in T law for variations of C in 
electron gas at low enough temperatures. However this approximation gives 
wrong results at T — > for systems with degenerated but partially filled (at 
T = 0) levels and at T — > oo for all systems. We do not apply this approxi- 
mation in our calculations and show wrong results for C given by this method 
only in some figures in subsequent sections. 

The polynomial method developed in our work, Ref. [T7], makes canonical 
{CE) calculations of C not more complicated than calculations corresponding 
to the EGCE description. Therefore in all figures given below we represent 
only the canonical results though in several cases the EGCE heat capacities 
are given to demonstrate to what extent C{CE) and C{EGCE) differ. 

Since the canonical occupation distribution n{e) is much closer to the step- wise 
distribution than the Fermi-Dirac one (see e.g. Refs. pT]|16|21] ) calculations 
of C with n{e) and f{e) at low temperatures lead to different but like results. 
As a rule EGCE calculations overestimate values of C at low T though with 
increasing T this difference disappears. However the simplicity of the basic 
relationship for C in the EGCE method gives the unconditional advantage to 
the EGCE method over the canonical one, Eqs. ([5]), ([6]), in interpretation of 
results. This feature of the EGCE approach will be exploited in subsequent 
sections. 

The amount of levels contributed markedly to the heat capacity increases with 
T giving rise to the growth of C. However the character of this growth is not 
monotonous and in the temperature evolution of C from zero temperature 
to values of order of or more than the Fermi energy (ep) four stages can be 
distinguished. Fig. [D 

• The low temperature (T < 6f) exponential growth of C up to a resonance 
temperature Tq ( a local maximum in C for some particle numbers). 

• A transition temperature region of inside which the rapid increase of C is 
replaced by practically linear T-dependence of C. 

• The linear ( or quasilinear ) growth of C with T which is well known for the 
electronic heat capacity in bulk metals. This stage is continuously trans- 
formed into the next stage. 

• The final stage of the evolution of C leading to the saturation: the heat 
capacity stops increasing with temperature and reaches the Boltzmann- 
Maxwell limit. 

To illustrate the characteristic features of C at each stage of the evolution 
we have performed canonical calculations of C for diverse systems mentioned 
in the Introduction in a wide range of N {N < 10^). In these calculations 
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as an energy unit for measuring single-particle energies and the temperature 
we employ the Fermi energy Sp of the considered system with a fixed A^. 
Therefore our data are appropriate for any fermion system with the same A^, 
shape, temperature and the type of confinement, however representation of 
these data on the absolute temperature scale requires the knowledge of the 
absolute values oi ep for each N. We believe that it is a detached problem 
and do not consider it here with the exception of a special model in Sec. 5 that 
gives ep{N) as a function of with respect to the corresponding bulk value, 
ep{bulk). 

Various types of confinement potentials we consider can be labeled by a quan- 
tity 7. For systems in cavities with reflecting walls (rectangular potentials) 
7 = -D/2, D being the spatial dimension, for all oscillator potentials 7 = 
This quantity determines the specific heat in the classic limit, at T/ep > 1 
{C/ksN = 7). For ^ 1 the value of 7 carries information concerning the 
type of confinement at linear variation of C with T {epN~^ T ep ) i.e. 
in the Sommerfeld heat capacity 

Cun/ksN = 7^—. (7) 

O Bp 



Sometimes, to stress peculiarities in C of a system with a finite A^ we adduce 
reduced values of the heat capacity C{N) /Cun- 

In some cases 5p, the mean level spacing for the Fermi gas with A^ ^ 1, is a 
more convenient unit to measure the temperature. : 

dp 

Sp = — — -, po{£p) = -fN/ep (8) 



Sp includes the level degeneration dp. For nonsymmetric systems dp = 2, for 
symmetric systems dp > 2. In particular, for spherical cavities the averaged 
degeneration dp is adopted to be equal to V^, Ref. [23], Po(^f) being the 
smooth part of the level density at ep and A^ 1 (more complicated depen- 
dence of po{ep) on A^ for small A^ discussed in Sec. 5 is not essential here). 

To simplify terminology we shall apply following names for systems under 
consideration. The term "system in a cavity" will imply a system of particles 
in a rectangular potential. If a cavity has a shape of parallelogram we shall 
distinguish two cube and a briquette. The latter is a parallelogram 

with different lengths of lateral ribs {L^ ^ Ly ^ L^). The inscriptions in figures 
"sphere" , "cube" , "briquette" and so on imply that the values of C are given 
for systems in cavities of corresponding shape, the inscription "oscillator" 
marks oscillator systems. 
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3D - oscillator N = 40 

2D - oscillator 




Fig. 1. The temperature dependence of the specific heat of 40 fermions in 3D-, and 
2D- isotropic oscillators (7 = 3; 2, respectively) and spherical and circular rectan- 
gular wells (7 = 3/2; 1, respectively). Four stages of the temperature evolution of 
C are distinguished. The first stage is the region of rapid temperature alterations of 
C/ksN, it includes a local maximum in C. The second stage is a transition region 
from the rapid increase to the linear (or quasilinear) growth with T that is the third 
stage. The fourth stage (not shown in the figure) is the region of the saturation and 
transition to the classic limit of C {C /ksN = ^). 
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3 Low temperature variations of the heat capacity 



For macroscopic fermion systems the first two stages in the evolution of C are 
absolutely inessential as the temperatures (T < 6p, 6f ~ £f/N) corresponding 
to these stages at ~ 10 are so low that practically inaccessible. However 
for mesoscopic systems the first stage is of the most interest since the variations 
of C v.s. T can give information concerning single particle levels near ep a.t 
temperatures accessible to measuring. (The absolute values of the temperature 
corresponding to the term "low temperatures" are determined by and the 
Fermi temperature which is material dependent: for metals ep ^ 10^ -r- 10^ K 
while for heterostructures ep ^ 10^ K and in trapped Fermi gases ep ^ IfiK.) 
Such possibility is opened owing to the properties of the function ip{e), Eq. IQ. 
These properties are determined by two factors constituting (p{e); the first one 
is the bell-wise function 



g/3(£-A) 



1 + e'^^^-") 



/(£)[!- /(e)] =T^, (9) 



A, the second factor is a parabola /3'^{e — X—/3^y with 



A + /3||. Thereby (f{e) is a two humped function with 



with a maximum at e 
a zero minimum at e 
the distance between two maxima ~ 5T, Fig. [2l Such form of (p{e) confirms 
Grimvall's conclusion, Ref. pi], that the heat capacity probes the level density 
in an interval ~ lOT. 



This distance weakly depend on /?|^. 



If > the left maximum is higher 



than the right one and vice versa at < 0. At high temperatures the left 
maximum does not take part in forming C as A <^ and /5A + is positive 
and approximately equal to 7 (see Sec. 6). 



In order to understand how ip{e) forms values of C at low temperatures it is 
useful to find the position of A with respect to the Fermi level (F). For this 
purpose we consider a model including only two levels: F and F + 1 supposing 
that Ep — Ep^i and ep+2 — ^f+i > £f+i — £f (the notations F ±k correspond 
to the k level above (+) and below (— ) F respectively). 

The position of A strongly depends on np, the occupation number of F at 
T = 0: 



1 < Hp < dp — 1] A 




Ep; 



(10) 



np = dp] A 



(Ep + Ep+i) +Tln 



d 



dp 



dp 



{Ep + Ep+i) /2.{n] 
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Fig. 2. The two humped function ip{£,t), Eq. Q, at f3d\/d(3 = and two temper- 
atures: t = 0.01 and t = 0.1 {t = T/ep). 
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Corrections to Eqs. fllOlllip are proportional to exp{—j3 \ eF±i — |) and 
omitted here. 

Eqs. f fTOfTTl) indicate that transition from the open shell to the closed one 
is accompanied by a drastic change of (3d\/d(3 determining the minimum . 
Whereas at np < dp {dp > 2) this quantity coincides with ep, i.e. the Fermi 
level does not contribute to C, at np = dp the value A + jSdX/dp is shifted 
to {ep + ep+i)/2, i.e. F-level takes part in forming C { in both cases pdX/dp 
does not practically depend on T). 

It implies that this model {EGCE) predicts an attenuating factor, exp[—0.5(3{ep+i — 
ep)], in the quotient of C{np) and C{dp): 



C{np)/C{di 



-0.5a; 



Id 



F+l 



2ni 



dp dp — np 
C{dp) = le 



np < dp — 1; 



0.5x^2 



dp+id} 



(12) 
(13) 



X = (3{ep+i - ep). 



The canonical approach to the same model does not give such temperature 
attenuation but confirms the strong dependence on np 

C{np- 0<np<dp)c^ ^!£^£±L^e-^x^ (14) 

dp — np + 1 



i.e. C{np = dp) prevails over C{np < df). Thus, C as a function of N has to 
oscillate with and reveals maxima at corresponding to closed shells. 

Heating extends ip{e) due to increasing the exponential factors in Eqs. ( |T3|T4l) 
and approaching other levels ( in the first place F — 1 and F + 2) to the 
maxima of (f{e) that increases their contribution to C. These two factors 
stimulate growth of C. 

A rather unexpected but quite natural consequence of the two humped char- 
acter of ip{e) is the appearance of a resonance against a background of the 
monotonous growth of C ( it is more evident in systems with high degenera- 
tions of levels). The resonance is conditioned by the coincidence of the points 
of the maxima in (p{e) with the energies of F and -F + l levels. Since 5T is the 
distance between the maxima of ip{e) the resonance temperature is determined 
by the difference ep+i — ep: 

To- ^(e^+i-e^). (15) 
5 
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(T/ep)N^ 



Fig. 3. Impact of the Fermi level filling on the heat capacity for shell = 40 in the 
spherical cavity, n^? is the particle number at T = on the Fermi level with degen- 
eration di? = 6. Results of the canonical (CE), equivalent grand canonical (EGCE) 
and grand canonical {GCE, A is temperature dependent, pdX/d(3 is omitted) meth- 
ods are displayed. It is obvious that the GCE method does not give correct values 
of C at T — > 0. In both other methods C — > at T — > 0. 



Such resonance amplification of C in the low temperature region was predicted 
as far back as the work of Frolich [1]. This local maximum can be found in 
the canonical calculations for ID- oscillator with ^ 1, Ref. [H], and for 
alkali clusters, Ref. p2j- However, the resonance nature of this phenomena is 
revealed here for the first time. 

Figs. [3], m [5] indicate that such maximum occurs in various systems . In all 
cases in these figures the temperature of the maximum practically corresponds 
to Eq. 1^. 
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Fig. 4. Top panel: The heat capacities of = 40 systems in 3-D-cavities of different 
shape (briquette, cube, cyhnder and sphere ). The dashed vertical lines mark values 
of t = (e_F+i — £f)/56f ( see Eq. (flSl) ) which practically coincide with the positions 
of the corresponding maxima in C. In this figure 5f = 'iep/'iN. Bottom panel: 
Fragments of the single-particle level schemes of sphere, cube, cylinder( the diameter 
is equaled to the height, the level degeneracy in cylinder is equal to 2 or 4) and 
briquette {Lx : Ly : Lz = I : l.l : 1.2, the level degeneration is equal to 2). 

In symmetric systems (spherical and cubic cavities, spherical oscillators) these 
temperature resonances reach maximum amplitudes for closed shell (i.e. in 
even systems) . However these resonances are high enough also ai np = dp — 
Development of such resonances with increasing np is shown for spheres and 
cubes in Figs. |3l [51 

If in the vicinity of F the spin degenerated {dp = 2) single-particle spectrum 
is uniform {Sp is the level spacing) then, as shown by Denton et al in Ref. [TT] . 
the local maximum in C at T ~ 0.25^ is explicitly expressed in even systems 
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Fig. 5. Variations of the reduced heat capacity (C/Cnn) vs the reduced temperature 
r = T/5p{cube) in the degenerate sheh of a cube, Spicuhe) = dp/ Po{£f)i Eq. ([8]). 
The averaged degeneration dp in a cube is adopted to be equal to 10 [dp = 10). 
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Fig. 6. Top panel: The reduced heat capacity (C/Cun) of A^-even and A^-odd systems 
vs T/6f for a briquette {L^ : Ly : Lz = 1 ■ 0.6e : tt), 6f = 4:£f/3N. Bottom panel: 
Fragment of the single-particle spectrum of the = 1440 briquette. 

while in odd ones at T ~ 0.45^? there is only irregularity in increasing C. 
However for realistic spectra level concentration near F can give maxima for 
both even and odd particle numbers. Fig. [6] indicates which variants can occur 
e.g. for briquettes with ~ 10^. Thus at a proper choice of A^ studying the 
temperature dependence of C can give the unique information concerning the 
difference ep+i — in mesoscopic systems. 

In systems with nonuniform spectra this difference can turn out to be small 
as compared with other level differences near F. In this case increasing T can 
give one more maximum in C. Examples of such cases are in Figs. [71 [8] which 
indicate that the second resonance is formed by not one pair of levels, as it 
takes place for the first resonance (STq = Aep+i^p = £f+i—£f) but in this case 
under the maxima of ip{e) there are at least two groups of levels (bunches). 
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Fig. 7. Top panel: The reduced heat capacity {C/Cnn) with two resonances for 
132 fermions in spherical cavity, ti and t2 are the positions of maxima. Bottom 
panel: The two humped functions ip{e,t), Eq. dH) at these temperatures. The first 
maximum is practically at ti = {sf+i —£f)/5£f but t2 corresponds to the presence 
of two groups of levels in the vicinity of the maxima in ip(e,t). Vertical lines mark 
the single-particle levels, their heights are proportional to the level degenerations. 



It should be mentioned that in the region of the first temperature resonance the 
divergence between C{CE) and C{EGCE) becomes maximal and then with 
increasing T this difference decreases and practically disappears at T > 2Sf, 
Fig.l 

To assess which differences AeF+i,F can occur in fermion mesoscopic systems 
and accordingly which temperature gives rise to resonances in C we have 
analyzed these differences in diverse systems with N ranged from 20 to 10^. 
These data are displayed in Fig. [in] in units Sf, Eq. ([H]). 
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Fig. 8. The same as in Fig. [7] but for the N = briquette. 

For the spherical oscillator the resonance temperature is strictly fixed (To ~ 
uj/2) but in spherical cavities there is a rather wide interval of possible val- 
ues of Tq. This interval would be much wider if one takes the same Sp (the 
denominator of x ) as in nonsymmetric systems i.e. without allowing for the 
averaged degeneration in spherical cavities. Possible additional degenerations 
have to be taken into account also at appearance of new types of symmetries 
arising when the ratios of the lateral lengths in briquettes or the frequencies 
in oscillator systems become equal to the ratios of small integers. 

For nonsymmetric systems in cavities the distributions of x = (£:_f+i —^f)/5f 
weakly vary with deformation of the system (excepting cases that lead to a new 
symmetry). However they evidently differ for 3D-cavities and 3D-oscillators. 
Nevertheless these distributions in both cases are rather wide. Therefore max- 
ima in C can occur at temperatures extending from very small values -C 6f 
up to (2 -h 3)(5i?. 
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Fig. 9. Difference between C{CE) and C{EGCE) as a function of T for briquettes 
{6f = ^ep/'iN) with different N. The most divergence is observed in the region of 
the resonance temperatures. 
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Fig. 10. Distributions (A) of values Aef+i^f/^f = i^F+i — £f)/Sf in different 
systems under consideration with 20 < < 10^. 6f is given by Eqs. ([8]). The 
histograms for A are built with the interval equal to 0.02. We use transcenden- 
tal numbers for the ratio of ribs in briquettes and oscillator frequencies to avoid 
accidental level degenerations. 



We would like to stress that the distributions in Fig. [TO] are not those consid- 
ered by Kubo and his followers. Kubo et al took into account the distributions 
of all spacings in an A^-particle system. Whereas in Fig. [10] the distributions of 
the only kind of spacings, {ep+i —£f) are presented but the particle numbers 
are in a very wide range 20<A^<10^at each type of confinement. 

After temperature resonances which complete the first stage of the evolution 
of C thermodynamic properties of a heated system are determined by several 
single-particle levels placed under maxima of ifiis), i.e. beginning with some 
temperatures the main role in forming C belongs to the smooth level density. 
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4 Transition to the regime of the smooth level density 



As mentioned above at T higher than the first or second resonance temperature 
several single particle levels (> 2) appear under and near the maxima of ip{e) 
that starts on averaging the level density p{e). This process is ended at the 
temperature Tgm, the temperature smoothing the level density oscillations, 
and accompanied by such interesting quantum size effect as oscillations of the 
heat capacity v.s. A^. This phenomenon results in the variations of the specific 
heat C/N v.s. A^, i.e. it exhibits one more paradoxical property of mesoscopic 
fermion systems in comparison with macroscopic ones. 

As well known the exact p{£), Eq. ([T]), can be represented as a sum of a smooth 
function of energy po{s) and an oscillating shell correction Sp{e) 

p(e)=po{e)+6p{e) (16) 



For simplest single-particle 3D-rectangular potentials values of poi^) were es- 
tablished many years ago and these results have been collected by Balian and 
Bloch [22j. For 3D-oscillators po{e) is given by Bohr and Mottelson [23]. The 
component Sp{e) is an infinite set of functions oscillating with e the period of 
which decrease as 1/n {n = 1,2, . . .). Each separate oscillation in 5p{e) with 
a period r„ will be named n-mode. 

In the previous section to interpret values of C we use the narrow distribution 
of ^{s) -function at low temperatures to choose some levels determining C. In 
this section because of dilation of ^{s) its role is changed and here it serves 
as an smoothing function to damp n-modes. In fact, if between the maxima 
of ip{e) there are ~ 3 periods r„, this mode will be suppressed and it will not 
practically contribute to C. It implies that the temperature T„ which removes 
n-mode and period r„ are connected by the condition 

5T„ ~ 3t:„, (17) 



i.e. Tn is of order of r„. If all r„ are proportional to a highest period tq, as it 
takes place in spherical oscillators 



P{e) 



ijj 



0J= 



i + 2E(-i: 



n=l 



27m 
cos e 



Tsmi smoothing all n-modes, is determined by this maximum period (for the 
spherical oscillator it is u;) 

Tsm ~ To. (19) 
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For asymmetric oscillators (all frequencies are different) the highest possible 
period amounts to the most frequency of cj^, ojy, u^- However if a system 
is moderately deformed {ur^ ujy ojz ^ uj = [ojxOJyUJzY/'^ the value of 
uo can serve as an estimate of such period. The oscillator frequency can be 
bound to physical parameters of the system by different ways. One of them 
consists in equating the oscillator Fermi energy with that for systems in a 
rectangular potential. In another way the empirical root-mean square radii (for 
3D-systems R ~ A^^/^, for 2D-systems R ~ A^^/^) are supposed to be equal to 
those calculated with oscillator functions ( [23] )• In both cases u ~ epN^^^'^ 
(SZ^-systems) or ~ epN"^^'^ (2D-systems). For closed oscillator shells with 
the main quantum number 'Rp and spin degenerated levels the exact relation 
between Kj^ and N is known: 3N = (i^F + 'i-){'^F + '2){'iiF + 3) for 3D-systems 
and = (KiT'-l- l)(Ki? + 2) for 2Z}-systems. Hence one obtains for 3> 1 that 
u ?a e^(3A^)"^/^ and EpiN)'^^'^ respectively. 

In spherical cavities the maximum period tq corresponds to oscillations cre- 
ated by supershells which are well studied j23f25f26ll?r] . Since the length of 
supershells ~ epN'^^^ period tq of density oscillations in mesoscopic spheres 
takes the same value. That is in agreement with the theoretical estimate of tq 



R being the radius of the sphere, R ~ N^^'^. 

Period tq for parallelepipeds is determined by the minimal lateral length of 
La;, Ly, Lz- For systems with ^ Ly ^ Lz ^ L = [L^LyLzY^^ the value of 
To is of the same order as for spheres: 

As far as high frequency modes are suppressed by heating the character of the 
A^-oscillations is altered. At small T, as pointed out in Sec. 3, C vs A^ at 
a fixed T demonstrates maxima at each A^ that corresponds to a filled level 
or shell. Amplitudes of these peaks are proportional to degenerations of adja- 
cent levels. The origin and temperature development of several such peaks are 
displayed in Fig. [TTJ The physical nature of the difference in peaks consists 
in variations of the Fermi level occupation numbers in adjacent systems at 
T 0. This difference has to be attenuated with growing temperature which 
stimulates regularization of the peaks due to removing high frequency modes 
in Sp{e). The gradual temperature refinement of the A^-oscillations is shown 
in Figs. [121 [131 [3 [El for diverse systems (spherical and deformed oscilla- 
tors, spherical and briquette cavities). These figures display such property of 



21 



(a) 



(b) 



1522 1760 

E shells M 

6 1760 

11 1734 (C) 

18 1688 

14 1614 

8 1556 

1 1522 



Fig. 11. Top panels: Temperature development of A^- oscillations of the heat capacity 
in spheres vs N, t = {T / ep)N^/'^ . Bottom panel: Fragment of the level scheme of 
the spherical cavity. N is the particle number of a system with the closed Fermi 
shell. L is the orbital momentum of the Fermi shell. 

mesoscopic systems (in comparison with macroscopic ones) as the dependence 
of C on the shape of the system ( even at a fixed type of confinement) and 
variations of the specific heat C/N with A^. 

As an illustration to the latter statement in Fig. [16] we compare the heat 
capacities of two systems containing = 2 ■ 10^ fermions, e.g. electrons in 
a metal grain. The first is a cube with the rib ~ (2 ■ 10^)^/^ro, tq ~ 0.1 
nm, the second system is a cube of the same size which is formed from 1000 
small cubes with 20 electrons in each. These small cubes are assumed to be in 
thermal contact i.e. they are kept at a common temperature but the mutual 
particle exchange is impossible. The top panel of Fig. [16] shows that at low 
temperatures the heat capacity of small systems, N < 10^, is practically 
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Fig. 12. N - oscillations of the canonical heat capacity C{CE) (bold line) and 
effective grand canonical one C{EGCE) (dotted line) at fixed T/u = 0.1, u is the 
spherical oscillator frequency. 
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Fig. 13. - oscillations of the reduced {C/Cun) heat capacity at different temper- 
atures r = {T/eF)N^/^ for 3D- deformed (top panel) and spherical (bottom panel) 
oscillators. Amplitudes of A^-oscillations in deformed oscillators are considerably 
lesser than in spherical ones. 

equal to zero and only behind iV ~ 10^ C gains the Sommerfeld (linear in 
T) values Cun- Therefore the canonical heat capacity corresponding to 10^ 
small cubes within the temperature interval < T/ep < 0.04 (the bottom 
panel of Fig. [T6ll remains considerably smaller than C of the "large" cube with 
= 2-10'^ electrons. Thus, regulating the size of mesoscopic grains in a sample 
gives a possibility to profoundly affect the thermodynamic characteristics of 
this sample at low temperatures. The bottom panel of Fig. [16] indicates once 
more that as a rule the EGCE method overestimates the values of C at low 
temperatures. 

In spherical cavities heating uncovers periodic structures oscillating with AA^ ~ 
N^/"^. These oscillations are closely connected with distribution of shells within 
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Fig. 14. N- oscillations of the heat capacity (C/Cun) in spheres at different reduced 
temperatures r = {T/£f)N^^^. 



25 



1 .2 




1 20 



Fig. 15. A^- oscillations of the heat capacity {C/Cun) in briquettes at different 
reduced temperatures r = {T /eF)N^/^ . 

supershells. Such bunches of shells can be occupied by large quantity of parti- 
cles ~ A^^/^, [23], accordingly the energy intervals between such level bunches 

As mentioned above (see Sec. 4) the values of C can be interpreted as the aver- 
aged values of the level density i.e. studying C v.s. one can get the averaged 
p as a function of N . Theoretically these quantities were calculated by means 
of averaging p{e) either with Lorentzians [27] or with Gaussians ^28j, in both 
cases the widths of these averaging functions were temperature independent 
and of order of CST^^. It is evident that at such temperature the result of 
averaging weakly depends on the type of averaging functions. 

In the frame of the EGCE formalism each oscillating functions (ri-mode) in 
p{s), i.e. cos(27r£:/rn + 0) for oscillators or cos{A.n ^/eps / Tk + 0) for cavities, 
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Fig. 16. Top panel: N- oscillations of the heat capacity (C/Cnn) in cubes. Bottom 
panel: The comparison of the heat capacity of a cube with N = 20000 with the 
total heat capacity of 1000 small cubes (each with = 20). 



(j) and (j) being independent phases, after averaging with (y9(£)-functions ( for 
the sake of simplicity we suppose here that A + pdX/dj3 ~ A) enters into C 
for oscillators or cavities respectively as cos (27rA/r„ + 0) or cos (47rA/r„ + 0). 
Thus the A^-oscillations of C in the EGCE approach (see Fig. [12]) is caused 
by variations of values of A/r„. For those temperatures, when the oscillations 
could be observable, A ^ that for large 3D-systems (A^ ^ 1) results in 
A/r„ ~ N^^^. Therefore in all figures the A^-oscillations in C for 3D-systems 
are presented as functions of N^/^. 

All arguments we have used for 3Z^-systems can be repeated for 2Z}-systems. 
The difference is only in the A^-dependence of oscillating periods tq and r„. 
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For 2D-systems 



T"0 ~ —r=] Tsm ~ To, (22) 



and the A^-oscillations have to be functions of N^^"^: These conclusions are 
confirmed by examples of the iV-oscillations of C in 2D isotropic oscillators ( 
Fig. [T71) and in 2D-systems with rectangular potential such as circles (Fig. [T8l) 
and rectangles (Fig. [T9|) . The suppershell pattern in circles is pronounced 
not so distinctly as in spheres and has a shorter beating period so that at 
r = (T/ep)N^^'^ > 0.5 the A^-oscillations in circles and rectangles are practi- 
cally identical. For circles the oscillating structures in the level densities are 
considered in Ref. [21] ■ 

This process which embraces beginnings of the A^-oscillations in C, appear- 
ance of periodic structures and their subsequent damping with T is in parallel 
to the de Haas-van Alphen effect in mesoscopic systems considered in our 
work, Ref. [2H], where variations of C and magnetic susceptibility in a linearly 
increasing magnetic field H were demonstrated at different temperatures. The 
analytical aspects of the analogy between C vs N and C vs were consid- 
ered in Ref. [30] . 

At heating beyond T^^ the heat capacity is determined only by po{s), the 
smooth level density, and the difference in occupation numbers does not affect 
the iV-dependence of C. This at least for very large systems, ^ 1, has to 
result in the T-linear regime of the heat capacity variations. In fact, see Fig.[20| 
for oscillators with A^ > lO'^ the reduced heat capacity C/Cun at T ~ Tsm is 
approximately equal to 1 whereas for lower values of A^ in oscillators and for 
all A^ including 10^ in cavities C/Cun < 1. 

For small A^ (A^ < 10^) temperatures of order of turn out to be so high, 
Tsm{N = 10^) ~ Bf/5 that the condition determining the T linear regime 
for C (T <C is not implemented. Practically after and near Tsm such few 
particle systems start going over to the classic limit of C (this transition regime 
we called in Sec. 2 as the quasilinear one). In systems confined in cavities with 
A^ > 10'^ the heat capacity in a temperature interval T > Tsm displays the 
T-linear dependence (Fig. [201) . However for a fixed shape of the system the 
factor of proportionality in C/N varies with the particle number that testifies 
to the more complicated A^-dependence of po(^) as compared with the case 
Ar> 1. 
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Fig. 17. A'"-oscillations of the heat capacity {C/Cun) in 2Z)-oscillators at different 
reduced temperatures r = {T/£f)N^/'^. 
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Fig. 18. N- oscillations of the heat capacity {C/Cun) vs N^/"^ in circles at three 
reduced temperatures r = {T /£f)N^/'^ . 



5 Deceleration of increasing c in temperature range Tgm <T < Sp 



On going out towards the hnear dependence on T the heat capacity slowly 
evolves to the classic Boltzmann-Maxwell hmit. Consequently there exists a 
temperature range (T > Tg^) within which the growth of C decelerates to 
attain to the saturation. To find the value of T giving visible deviations from 
the linear low N 1 the smooth level density pQ{e) can be represented in 
the form: 



N 



,7-1 



(23) 
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Fig. 19. N- oscillations of the heat capacity {C/Cun) vs N'^/'^ in rectangles 
[Lx : Ly = 1 : e/vr) at three reduced temperatures r = {T / £f)N'^/'^ . 

then the integral for C, Eq. ([2]), with po(^) instead of Pex^^) can be calculated 
by the standard way up to terms (/Jei?)"^ and (/Ssi?)^ exp(— /Jei?) 

rih \^ (7-l)(9-27K 3(7-2)(7-3) , 



For few particle systems Eq. fl2^ can be used to obtain rough estimates of 
such temperature. The origin of the last term in Eq. flMl) is the finiteness 
of the value of ep/T that is usually supposed to be infinitely large. However 
it is not the case if {l3eF)~^ is of order of several tenths. Eq. fl2^ shows that 
(20-i-30)% deviations from linearity are observed at (Pep)'^ ~ 0.3 for 7 = 3/2 
( 3D-cavities) and at (Pep)'^ ~ 0.2 for 7 = 3 (3D-oscillators). Absolute values 
of this temperature corresponding to these deviations are material dependent. 
So for metal clusters it is ~ 10'^ -^ 10^ i^T while for heterostructures ~ 100 K. 
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Fig. 20. The reduced (C/Cun) heat capacity vs r = {T/£p)N^/^ (in the temperature 
region where C is T linear for large A^) for 3D- spherical (a) and deformed (b) 
oscillators and for spherical (c) and briquette (d) systems in rectangular cavities at 
different A^. 
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Fig. [2T] displaying the specific heat in cavities demonstrates the dependence of 
the specific heat on A^. This is the obvious manifestation of the comphcated 
iV-dependence of the smooth level density po{e) that was mentioned in the 
end of the previous section. 

For 3D-cavities, 7 = 3/2, po(^) consists of the volume {V), surface (S) and 
linear terms [22] 

V (h^Y''^ . s fh^ Y' (h^ y 1 



The explicit dependence of po{e) on the sizes of the system can be avoided 
due to the relationship between ep and a parameter of length inherent in 
the system, Eq. fl2B]) . Such parameter for spheres is the radius R, for cubes 
the lateral length L. The same quantity L is suitable for briquette if L = 



where Xp for spheres is a Bessel function root, for briquettes 



X| = tt' (nl/al + nl/al + nl/a'^ , (27) 
ax : ay : Uz = : Ly : Lz, a^ayaz = 1, 

For cubes a^ = ay = a^ = I. 

Then, the replacement of /2mR^ or ^^/2mL^ by ep/Xp reduces Eq. (1251) 
to 

Po(^) = — E 7S — • (28) 



Coefficients straightforwardly found by comparison of Eqs. (!25|) . (126!) 

and Eq. 



~ Xp/N. (29) 



Calculations of C with po(^); Eq. (128|) . lead to the result similar to Eq. (^^, 
but the total factor determining the linearity in T and factors before (/Se^?)"^ 
and (Pep)'^ exp{Pep) gain a more complicated dependence on as compared 
with Eq. (IMl) . In particular instead of 7 = 3/2 (the first factor in Eq. (IMl) ) 
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now we have I]^7a^(^^j . (The quantity Aq is defined by the equation 

^7^7(7^)^ = 1). For small (< 10^) this factor > 3/2, that is caused 
mainly by coefficients 03/2, Eq. fl29l) . which for small is more than 1. For 
A^ ^ 1 the dependence Xp on A^ can be established by using the integral 
indep ending of 

1 

— / p{e)de = Y^a^. (30) 

T 



For large values of A^ Eq. fl30l) is practically equal to 1. The surface and linear 
term in p{e) aX N > 10^ are unessential and so from Eq. fl5Ul) it follows 



= 9ttN/4: (sphere) or Sn'^N^cube or briquette). (31) 



However at small A^ the negative surface term (~ Xp) plays the important 
role, raising the values of Xp as compared with Eq. fl3Tl) . This explains exceed- 
ing oi C /N for small A^ over large ones at the same temperature in Figs. 



The values oi C/N for small N can be affected even by the linear term in 
Po{s). The surface and linear terms in C have different dependence on shape 
parameters of briquettes, Eq. fl27|) : the surface term ~ Z^Q^i"^ and linear term 
~ These sums can be essentially distinguished for strongly stretched 

and flattened shapes. That can give rise to difference in C at small A^. 

Fig. [22I where C/N for a spherical oscillator are displayed, shows that A^- 
dependence of C/N for oscillators at T ~ e^? are practically absent since for 
this type of conflnement, 7 = 3, the additional term in p{e) is much less than 
the flrst (~ Ar-2/3) 



Po{e) = {u^UyUj^) 



-1 



2 / 2 I 2 I 2\ 



N-^/'\ (32) 



Only for small A^ and strongly deformed oscillators the second term in Eq. fl32|) 
could play some role. 

Thus, the variations C/N ysT aiT > Tg^ indicate that the factor determining 
the linearity in T can get the A^-dependence caused by the type of conflnement. 
Probably realistic potentials are something middle between rectangular and 
oscillator potentials. Therefore studying the speciflc heat could give additional 
information concerning the type of conflnement in a fermion system under 
consideration. 

Variations of e^? vs A^ in fermion mesoscopic systems could be predicted by a 
theory ascertaining connection of A^ with linear sizes of the system. Here we 
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will consider the jellium model postulating that the linear sizes are propor- 
tional to 

R = roN^/^ [sphere], L = loN^^^briquette), (33) 



To, Iq being independent of A^. 

Then, applying Eqs. fl26l) . fl3T|) . fl33l) one obtains 



(34) 



^(bulk) JY2/3 \^g7]-y 
£i.(Ar) X2 / 1 \2/3 



(briquette), (35) 



where Ep^'"''^ is the bulk Fermi energy. 

The results of calculations of EpiN), Eq. (|34l) . for spherical systems are dis- 
played in Fig. [251 In the same figure the smooth curve is plotted. It is obtained 
by replacing the exact values of Xp by approximate ones found by using 
Eq. fl30|) (assuming that J^o^-y = !)• This curve allowing for only the surface 
term (see also |32j) in pq, Eq. fl25|l . gives the following smooth dependence 

(smooth) 7> T 

Ep on N: 



(bulk) 

c 

{smooth) 



{smooth) 

(36) 



ep 



Ep 



{bulk) \3N J 




(^) (E«~')'^' (briquette), (37) 



the definition of is given by Eq. (l29ll . Fig. [25] and Eqs. (1361) . (137]) show that 
in the model in which the parameter of spatial density (r^ or Iq), Eq. (1551) . is 
independent of iV the energy oscillates with N around 5^^'"°°*'') and even 
for small particle numbers (N ~ 50) e^^^ deviates from e^"''^^ not more than 
by 40%. 



6 The classic limit of C 



Though in the temperature range T^^ < T < Ep the specific heat (C/N) 
reveals some dependence on N, in the limit of superhigh temperature T/ep ^ 
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Fig. 21. The specific heat in a cube (7 = 1.5) for different N at high temperatures. 
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Fig. 22. The specific heat in a spherical oscillator (7 = 3) for = 40 and 1000 at 
high temperatures. 

1 this dependence completely dies out and irrespectively of the shape and 
particle numbers C/N trends to the classic Boltzmann-Maxwell limit 

C/keN = 7 (38) 



(remind that 7 = D/2 for cavities and 7 = D for oscillators). The comparison 
of C obtained by three methods {CE, GCEaxid EGCE, see Sec. 2) in the high 
temperature region up to is given in Figs. [231 EH It is obviously that at 
high temperatures and at T — > 00 for all systems EGCE gives practically 
the same results as CE and both methods lead to the correct classic limit 
while GCE gives wrong results for C at T > Q.2ep. Consider the behavior of 
the specific heat in the limit of superhigh temperature, T/sp ^ 1, in more 
details by using the EGCE formalism and begin with the calculation of to 
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Fig. 23. Comparison of the values of C calculated by three methods: CE^ EGCE 
and GCE (see Sec. 2). Calculations are performed for systems with 40 fermions 
confined in rectangular potentials of different spatial dimensions. On the scale of 
the figure C{EGCE) coincides with G{CE) at all temperatures. At T > {).2eF the 
GCE method overestimates the values of C and gives wrong results at T — > oo. 

find pd\/dp 



CXJ 



(39) 



The simple form of Eq. (139!) takes place since at T ^ e^? the chemical potential 
becomes negative A < ( the energy scale is chosen so that all single particle 
level energies are positive). Taking advantage of the expansion of po(^) in 
powers of e one can transform Eq. (I39p into the following 
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Fig. 24. The same as in Fig. [23] but for the isotropic oscillators. On the scale of the 
figure C{EGCE) coincides with C{CE) at all temperatures. 



(40) 



Eq. (HOl) indicates that the more is the temperature the less terms with 7 < 
Imax as compared with the first term with 7 = 7maa;- Below index "max" will 
be omitted. In this approximation Eq. (HUI) defines [3\ and [3'^dX/d[3: 



/3X 



T 

7 In h In + In r(7 + 1) 



A + /3^1-e^V.r(7 + !)(-) ^7. 



(41) 
(42) 
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Fig. 25. The Fermi energy ep (solid line) in spheres vs N in the jellium model, 
Eq. (j34p . where tq is material dependent parameter independent of A^, {R = tqN^/^, 
ro is the radius of a sphere containing one particle), g^^™""*^) ig tj^e smooth Fermi 



energy, Eq. (j36j) . 



Comparison of Eq. (jUj) and Eq. (H2l) emphasizes the importance of the term 
dX/dp because only taking into account this term makes the quantity /5(A + 
PdX/dP) independent of T. 

Also as in Eq. ( HOl) the calculations of C is performed with the only term in 
p{e) with maximum 7. 



C/kB = N-^a^e^^Ep^ J{x- -ffe^-^e-'' 



de 



(43) 
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In this equation x 
Eq. (go]). 



(3e and ^(A + (3§] 



is replaced by 7 in accordance with 



N-fa,e''\(3eF)-^T{^) [(7 + 1)7 



27" + 72 



(44) 



Comparison of first factors ( before square brackets) with Eq. fHOj) gives the 
final result, Eq. (138|) . The same result can be obtain in the canonical formalism 
starting with the canonical partition function Zj\f. As shown in Ref. [T7j Zjv 
is equal to [A^], a symmetrical polynomial of the power N in variables = 
exp(— /Se^), Es being the single-particle energy of state s. According to the 
known Warring formulae each symmetrical polynomial can be represented 
through power sums 5*^ 

00 

Sn = T, dsq: = I Pe.{e)e''^'de, (45) 



at high temperatures Pexi^) — ^ Po(^)- We write out only two first terms of 
this representation to show that the second term and all others disappear at 
T — ^00 



For superhigh temperatures: 



Sn = -s(n/?£^)-^r(7), (47) 
n 

S2/{S^f = (/?£f )VAr2^+ir(7) It-^oo— 

Thus, to find C in the high temperature limit it is sufficient to take only the 
first term in Eq. P6|) into account and use Eq. (H7|) at n = 1 that straightfor- 
wardly gives rise to the final result, Eq. fl38l) : 

C/kB = (3'^ \n[N] = 13'-^ {N\nS,} = ^N, (48) 



since In Si = —7 ln.(3 + terms independent of /3. 
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7 Conclusion 



The investigation of the heat capacity variations at heating from T = up to 
T Sp allows us to infer what information can be obtained at measuring the 
fermion mesoscopic heat capacity. 

At temperatures of order of or lesser than the average level spacing 5 p (i.e. 
dpSp/lN where dp is the average level degeneration near the Fermi level) 
experiments make it possible, as it arises from Section 3 and 4, to establish 
single -particle level structure in the vicinity oisp. For this purpose two kinds 
of experiments could be performed on systems of the same type with a small 
variations in particle numbers. The first one is studying low temperature lo- 
cal maxima or irregularities in the heat capacity temperature increase that 
gives the resonance temperature directly connected with Spj^i — £p (Sect. 3) 
in each of such systems. The second kind of experiments on the same systems 
is measuring C as a function of N either at a fixed temperature smaller than 
5p or at temperatures of the same order but decreasing with N (~ N~^^^ for 
3D-systems or ~ iV"i/2 for 2D -ones). These experiments could give maxima 
in C for those values of that correspond to filled levels or shells. These ex- 
perimental results (C v.s. T for adjacent N or C v.s. N around some value Nq) 
can be used to reconstruct the energy level sequence near ep corresponding to 

As discussed in Sections 2 and 4 the heat capacity can be interpreted as the 
temperature averaged level density. Studying C v.s. N at T 5p gives too 
much details in the level density as a function of A". However at T ~ 0.57^^ 
( behind Tsm all level density oscillations discontinue, Tgm — epN~^/^ for 3D- 
systems and Tgm — SpN^^^"^ for 2D-oncs) measuring C uncovers the main 
periods of oscillations in C and accordingly in the level density that can be 
applied to researching other mesoscopic phenomena depending on the level 
density. 

Besides the theoretical interest studying C v.s. A" can find practical appli- 
cations. For example, the heat conductivity is known to be proportional to 
the heat capacity and so variations of C have to be taken into account if a 
nanogranule material is chosen as a thermoinsulator at very low temperatures 
when the electron heat capacity prevails over the lattice one. At such tem- 
peratures, as shown in Sec.4, the electron heat capacity of some mesoscopic 
systems, i.e. at some values of the particle number, takes practically zero values 
thereby such nanogranule materials are also devoid of the heat conductivity. 

This work is supported by the ISTC under grant Nr. 3492. The authors are 
much indebted to V.P. Chechev and R.B. Panin for the help in the work. 
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